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The ring closure reaction of C 4 H 3 with acetylene to give phenyl radical is one proposed mechanism 
for the formation of the first aromatic ring in hydrocarbon combustion. There are two low-lying 
isomers of C 4 H 3 ; l-dehydro-buta-l-ene-3-yne (n-C 4 Hj) and 2-dehydro-buta-l-ene-3-yne 
(iso-C 4 H 3 ). It has been proposed that only «-C 4 H 3 reacts with acetylene to give phenyl radical, and 
since iso-C 4 H 3 is more stable than / 1 -C 4 H 3 , formation of phenyl radical by this mechanism is 
unlikely. We report restricted Hartree-Fock (RHF) plus singles and doubles configuration 
interaction calculations with a Davidson s correction (RlfFT I -t-2-l-l2) using the E^unning 
correlation consistent polarized valence double zeta basis set (cc-pVDZ) for stationary point 
structures along the reaction pathway for the reactions of /1-C4H3 and iso-C 4 H 3 with acetylene. 
"-C 4 H 3 plus acetylene (9.4) has a small entrance channel barrier (17.7) (all energetics in parentheses 
are in kcal/mol with respect to iso-C 4 H 3 plus acetylene) and the subsequent closure steps leading to 
phenyl radical (-91.9) are downhill with respect to the entrance channel barrier. ISO-C4H3 plus 
acetylene also has an entrance channel barrier (14.9) and there is a downhill pathway to 
1-dehydro-fulvene (—55.0). 1-dehydro-fulvene can rearrange to 6 -dehydro-fiilvene (—60.3) by a 
1,3-hydrogen shift over a barrier (4.0), which is still below the entrance channel barrier, from which 
rearrangement to phenyl radical can occur by a downhill pathway. Thus, both /1-C4H3 and iso-C 4 H 3 
can react with acetylene to give phenyl radical with small barriers. © 1995 American Institute of 
Physics. 




I. INTRODUCTION 

The formation of the first aromatic ring (phenyl radical) 
is widely believed to be the rate limiting step in the forma- 
tion of soot in hydrocarbon combustion. One mechanism for 
the formation of phenyl radical involves the stepwise addi- 
tion of acetylenes.^ One possible pathway is the reaction of 
the vinyl radical with acetylene to give C 4 H 3 followed by 
subsequent addition of another acetylene and ring closure to 
give phenyl radical. Another pathway to C 4 H 3 is the reaction 
of vinylidene with acetylene to give vinylacetylene 
(buta-l-ene-3-yne),^ which can be converted to C 4 H 3 (1- or 
2-dehydro-buta- 1 -ene-3-yne) by abstraction of a hydrogen. 
Note that reaction of acetylene with itself does not occur by 
any low energy pathway, and vinylidene is 43 kcal/mol 
above acetylene;^’^ thus, the formation of vinylacetylene 
from two acetylenes requires a moderate activation energy. 
The structure and properties of C 4 H 3 have been studied by 
Ha and Gey"^ using PUMP4/6-31G** theory. They find two 
low-lying isomers l-dehydro-buta-l-ene-3-yne (/ 2 -C 4 H 3 ) and 
2-dehydro-buta- l-ene-3-yne (iso-C 4 H 3 ). ISO-C 4 H 3 is found to 
be 7.3 kcal/mol below «-C 4 H 3 (including zero-point effects). 
Miller and Melius^ also noted that iso-C 4 H 3 is thermody- 
namically more stable than and thus most of the 

C 4 H 3 should be in the form of the iso isomer even at com- 
bustion temperatures. Thus, it is of interest to study the re- 
action of both iso-C 4 H 3 and / 1 -C 4 H 3 with acetylene, as dis- 
cussed herein. 

In Sec. II the technical details of the calculations are 
discussed, while Sec. Ill discusses the results, and Sec. IV 
concludes the paper. 


II. COMPUTATIONAL DETAILS 

Two different basis sets were used in these calculations. 
The stationary points were located with restricted Hartree- 
Fock (RHF) derivative calculations using the valence double 
zeta set of Dunning and Hay.^ The basis set for C is the 
{9s5p)f[2is2p] basis and the H basis is i4s)/[2s], i.e., the 
polarization functions are omitted. The subsequent internally 
contracted configuration interaction (ICCI) calculations were 
based on the RHF reference wave function and used the 
Dunning correlation consistent polarized valence double zeta 
basis set.^ 

The geomedy optimizations used the sirius/abacus 
system of programs,® while the ICCI calculations were car- 
ried out with MOLPRO.^’*^ All electrons were correlated ex- 
cept for the C 1 5 like electrons. A multireference analog of 
the Davidson’s correction*^ was added to the ICCI energies 
and is denoted by +0. 

III. DISCUSSION 

Table 1(a) shows the computed ICCI energies of the sta- 
tionary point structures for the reaction of iso-C 4 H 3 with 
acetylene, while Table 1(b) shows the same information for 
the reaction of n-C 4 H 3 with acetylene. The computed har- 
monic frequencies are given in Tables Ala and Alb of the 
PAPS material.*^ The zero point effects were estimated as 
1/2 the sum of the harmonic frequencies and these are in- 
cluded in the relative energies, which are given in the last 
column of Tables 1(a) and 1 (b). Finally, the Cartesian coordi- 
nates for the stationary points are given in Tables Alla and 
Allb of the PAPS material. 
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TABLE I. (a) Computed energies for stationary points on the iso-C 4 H, plus 
C 2 H 2 surface."* (b) Computed energies for stationary points on the n-C 4 H^ 
plus C 2 H 2 surface.*’ 


Geometry 

ICCI (ICCI + e+230.) 

ZPE 

A£ 

reac 

(a) 

-230.551 40(- 0.685 46) 

0.075 38 

0.0 

spl 

-230.530 32(- 0.668 48) 

0.082 08 

14.9 

mini 

-230.612 79(-0.742 74) 

0.087 99 

-28.0 

sp3 

-230.602 l4(-0.733 03) 

0.085 70 

-23-4 

min3 

-230.612 50(- 0.742 39) 

0.087 84 

-27.9 

sp6 

-230.586 35(-0.719 08) 

0.086 60 

-14.1 

min6 

-230.657 85(-0.789 94) 

0.092 23 

-55.0 

sp7 

-230.553 I5(-0.688 29) 

0.084 52 

4.0 

min 7 

-230.665 83(- 0.797 63) 

0.091 41 

-60.3 

sp8 

-230.601 10(-0.734 03) 

0.086 29 

-23.6 

min8 

-230.610 98(-0.741 16) 

0.087 99 

-27.0 

sp9 

-230.529 84(-0.664 82) 

0.082 95 

17.7 

min9 

-230.659 85(-0.791 56) 

0.091 76 

-56.3 

splO 

-230.587 26(-0.720 29) 

0.086 40 

-14.9 

minlO 

-230.618 72(- 0.748 86) 

0.086 78 

-32.6 

reactants 

(b) 

-230.537 98(-0.671 01) 

0.075 87 

9.4 

sp2 

-230.530 09(- 0.663 57) 

0.081 72 

17.7 

min2 

-230.605 13(- 0.737 60) 

0.088 29 

-24.6 

sp4 

-230.595 01(-0. 728 48) 

0.085 96 

-20.4 

min4 

-230.610 98(-0.741 16) 

0.087 99 

-27.0 

sp5 

-230.589 19(-0.729 47) 

0.087 62 

-19.9 

min5 

-230.720 30(-0.850 62) 

0.094 09 

-91.9 


“Energies in a.u, except for A£ which is the relative energy in kcal/mol with 
respect to iso-C 4 H 3 +C 2 H 2 including zero point energy. 

’’Energiees in a.u. except for AE which is the relative energy in kcal/mol 
with respect to iso-C 4 H 3 +C 2 H 2 including zero point energy. 


The stationary points structures for the reaction of 
iso-C 4 H 3 with acetylene are shown in Fig. 1, while those for 
the reaction of with acetylene are shown in Fig. 2. 

Figure 3 shows the stationary point structures for the conver- 


•Pl mini *p3 _ min3 


sp6 minS sp5 minS 


FIG. 1. Stationary point structures for the reaction of 
iso-C 4 H 3 +C 2 H 2 “>C 6 H 5 (phenyl). The following nomenclature is used. Mini 
and min3 = l-dehydro-3-methylene-penta- l-ene-4-yne. Min6= 1 -dehydro- 
fulvene. Min7=6~dehydro-fulvene. Min8=l-dehydro-hexa-L3-diene-5-ync 
=Min4. Min5-phenyl radical. Structures designated by min are minima on 
the PES while structures designated by sp are saddle points. This convention 
is also used for Figs, 2 and 3. 


sp2 min2 sp4 



min4 sp5 minS 

c^:) 

4 ^^ 


FIG. 2. Stationary point structures for the reaction of «-C 4 H 3 + C 2 H 2 -^C 6 H 5 
(phenyl). The following nomenclature is used. Min2 and min4= I-dehydro- 
hexa-l,3-diene-5-yne. MinS = phenyl radical. 


sion of min 6 to min 10. (Here structures denoted by min are 
minima and structures denoted by sp are saddle points on the 
potential energy surface.) The energetics for these processes 
are shown in Figs. 4-6, respectively. 

Consider first the entrance channel barriers for the reac- 
tion of r?-C 4 H 3 and iso-C 4 H 3 with acetylene. A 1 -C 4 H 3 plus 
acetylene is 9.4 kcal/mol above iso-C 4 H 3 plus acetylene, but 
the barrier to addition is 14.9 kcal/mol for iso-C 4 H 3 plus 
acetylene and 8.3 kcal/mol for / 1 -C 4 H 3 plus acetylene. The 
higher bamer for iso-C 4 H 3 plus acetylene probably results 
from larger nonbonding interactions between the two acety- 
lenic groups, which come closer together in the iso-C 4 H 3 
plus acetylene case. This is reflected in the larger angle be- 
tween the acetylenic groups in Fig. 1 as compared to Fig. 2. 

Addition of A 1 -C 4 H 3 to acetylene leads to min 2 . Inversion 
about the CH group of Cl interconverts min2 and min4 with 
a bamer of 4.2 kcal/mol. Min4 closes to min5 (phenyl radi- 
cal) with only a small barrier. Thus, this reaction path leads 
to phenyl radical with no bamer other than the entrance 
channel addition barrier. 

Addition of iso-C 4 H 3 to acetylene leads to mini. Mini is 
connected to min3 by sp3. The conversion from mini to 


sp9 niln9 

spIO minlO 


FIG. 3. Stationary point structures for the conversion of min6 to min 10. The 
following nomenclature is used. Min9=2-dehydro-fulvene. Min 10=2- 
dehydro-hexa- 1 ,3-diene-5-yne. 
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sp7 — > minlO 


40.0, 


sp9 


0.0 


-40.0_ 


' 80 . 0 . 



FIG. 4. Energetics from SDCI + 0 calculations with the cc-pVDZ basis set 
for iso-C 4 H 3 +C 2 H 2 — (phenyl). 


min3 involves an inversion of the CH group about Cl. This 
involves only a 4.6 kcal/mol barrier. Sp6 connects min3 to 
min6 (1-dehydro-fulvene). Sp7 is the saddle point for a 1,3- 
hydrogen shift which connects min6 (1-dehydro-fulvene) to 
min? (5-dehydro-fuIvene). This process involves a large bar- 
rier; sp7 is 4 kcal/mol above iso-C 4 H 3 plus acetylene, but 
still below the entrance channel saddle point, so this channel 
is energetically accessible. Min7 can convert to min8 via 
sp8. Min8 is the same structure as min4 in the «-C 4 H 3 plus 
acetylene reaction pathway. This can then convert to phenyl 
radical via a downhill pathway. Thus, both iso-C 4 H 3 plus 
acetylene and rr-C 4 H 3 plus acetylene can give phenyl radical 
by pathways which involve no barriers, which are larger than 
the entrance channel addition barriers. 

Sp9 is a saddle point connecting min? to min9 (2- 
dehydro-fulvene). Sp9 is 17.7 kcal/mol above iso-C 4 H 3 plus 
acetylene, which is still below the entrance channel saddle 
point. However, sp9 is a second order saddle point, i.e., it has 
two negative eigenvalues [see Table 1(a)] one of which leads 
to symmetry breaking. Thus, there may be a lower energy 
path to min9 if the symmetry is lowered. This point was not 
examined further in this work. SplO connects min9 to min 10. 


n-C4H3 + C2H2 — > C6H5 

40.0 

sp2 



- 120.0 


FIG. 5. Energetics from SDCH -0 calculations with the cc-pVDZ basis set 
for / 7 -C 4 H 3 +C 2 H 2 -+C 6 H 5 (phenyl). 


FIG. 6 . Energetics from SDCI+0 calculations with the cc-pVDZ basis set 
for the conversion of min 6 to min 10 . 

IV. CONCLUSIONS 

The ring closure reactions of l-dehydro-buta-l-ene-3- 
yne (n-C 4 H 3 ) and 2-dehydro-buta-l-ene-3-yne (iso-C 4 H 3 ) 
with acetylene have been studied using RHF plus singles and 
doubles Cl plus a Davidson’s correction with the Dunning 
cc-pVDZ basis set at stationary point structures determined 
by RHF derivative calculations using a polarized valence 
double zeta basis set. Zero-point effects were estimated 
based on the RHF harmonic frequencies. 

n-C 4 H 3 plus acetylene (9.4) has a small entrance channel 
barrier (17.7) (all energetics in parentheses are in kcal/mol 
with respect to iso-C 4 H 3 plus acetylene) and the subsequent 
closure steps leading to phenyl radical (-91.9) are downhill 
with respect to the entrance channel barrier. ISO-C4H3 plus 
acetylene also has an entrance channel barrier (14.9) and 
there is a downhill pathway to 1-dehydro-fulvene (-55.0). 
1-dehydro-fulvene can rearrange 6-dehydrofulvene (-60.3) 
by a 1,3-hydrogen shift over a barrier (4.0), which is still 
below the entrance channel barrier. 6-dehydro-fulvene can 
rearrange to phenyl radical via l-dehydro-hexa-l,3-diene-5- 
yne (-27.0). Thus, both «-C 4 H 3 and iso-C 4 H 3 can react with 
acetylene to give phenyl radical with small barriers. 

With respect to iso-C 4 H 3 plus acetylene, all of the open 
structures (mini, min3, min8=min4, and minlO) are from 
27.0 to 32.6 kcal/mol lower, while the dehydro-fulvene struc- 
tures (min6, min7, and min9) are from 55.0 to 60.3 kcal/mol 
lower, and the phenyl radical is 91.9 kcal/mol lower, respec- 
tively. Thus, since all these structures can interconvert with 
barriers below the entrance channel barrier, the ultimate 
product is expected to be phenyl radical, although other 
species may still play a role as reactive intermediates. 

ACKNOWLEDGMENTS 

S.RW. was supported by NASA Cooperative Agreement 
Number NCC2-478 to FLORET Institute and by NASA 
Contract No. NAS2- 14031 to FLORET. 

I. Classman, in Proceedings of the 22nd Symposium (International) on 
Combustion (The Combustion Institute, Pittsburgh, 1988). 

“S. P. Walch and P. R. Taylor, J. Chem. Phys. (in press). 

^M. M. Gallo, T. P. Hamilton, and H. F. Schaefer III, J. Am. Chem. Soc 
112, 8714 (1990). 

'‘T-K Ha and E. Gey, J. Mol. Struct. 306, 197 (1994). 


J. Chem. Phys., Vol. 103, No. 19, 15 November 1995 



8547 


Stephen R Welch: MEP for reactions of C 4 H 3 with acetylene 


A. Miller and C. F. Mellius, Combust. Flame 91, 21 (1992). 

T. H. Dunning, Jr. and P. J. Hay, in Methods of Electronic Structure 
^Theory, edited by H. F. Schaefer III (Plenum, New York, 1977). 

T. H. Dunning, Jr., J. Chem. Phys. 90 , 1007 (1989). 

® SIRIUS is a MCSCF program written by H. J. Jensen and H. Agren and 
ABACUS is a MCSCF derivatives program written by T. Helgaker, H. J. 
Jensen, P. J^renson, J. Olsen, and P. R. Taylor. 

^’H.-J. Werner and P. J. Knowles, J. Chem. Phys. 89, 5803 (1988). 

P. J. Knowles and H.-J. Werner, Chem. Phys. Lett. 145 , 514 (1988), 


" S. R. Langhoff and E. R. Davidson, Int. J. Quantum Chem. 8, 61 (1974). 
See AIP document No. PAPS JCPSA6- 103-8544- 18 for 18 pages of tables 
of supplemental material. Order by PAPS number and journal reference 
from American Institute of Physics, Physics Auxilary Publication Service, 
Carolyn Gehlbach, 500 Sunnyside Boulevard, Woodbury, New York 
11797-2999. Fax: 516-576-2223; e-mail: Janis@aip.org. The price is $1.50 
for each microfiche (98 pages) or $5.00 for photocopies of up to 30 pages, 
and $0.15 for each additional page over 30 pages. Airmail additional! 
Make checks payable to the American Institute of Physics. 


J. Chem. Phys.. Vol. 103, No. 19, 15 November 1995 



